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Abstract 

Based on quantum theory, we investigate the distribution of the electrons scattered in nonlinear 
Compton effect by an electromagnetic plane wave. The monochromatic case, examined in detail, 
reveals features of the electron distribution, useful in the understanding of the pulsed plane wave 
case for particular intensity and electron energy regimes. The graphs displayed focus on the case 
of head-on or near head-on collision of an energetic electron with an electromagnetic circularly 
polarized pulsed plane wave and show that the deviation in direction is extremely small, while the 
distribution in energy can be visibly different from that of the initial electron. Two pulse shapes, 
several laser intensities and high incident electron energies are considered. 
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I. INTRODUCTION 



Nonlinear Compton scattering (NLCS) is one of the simplest processes predicted by 
quantum theory [IH3] and was detected in the head-on collision of an energetic electron 
beam with an intense laser beam [H |5]. 

In the case of a pulsed electromagnetic plane wave (a plane wave with a finite extension in 
the direction of propagation), the model we use for the laser beam, an electron distribution 
at the end of the pulse different from the initial one can not be predicted within classical 
electrodynamics (CED) without including the radiation reaction (RR) as, according to this 
theory, the pulse leaves each electron with the same momentum it had at the beginning of 
the pulse. On the contrary, the emission of radiation can be described by CED as Thomson 
scattering: the charged particle accelerated by the electromagnetic field emits radiation 
during its well determined motion. The inclusion of RR in CED takes into account the 
energy lost by the electron and leads to a final electron momentum different of the initial 
one. The classical description of RR was analyzed in several publications in which not only 
the formalism was discussed but also the effects of radiation reaction on Thomson scattering 
(references can be found in the very recent review [6]). In quantum theory RR was considered 
only very recently [Z],[8]. It is argued that the standard description of NLCS, as a single one- 
photon emission, using the Volkov solutions of Dirac equation, does not include radiation 
reaction and that the mechanism for it is the incoherent multiple one-photon emission by 
the electron. 

Up to now comparisons between quantum and classical predictions were done only for 
the emitted radiation spectrum. In the work presented in this paper we do no not include 
RR effects. 

In quantum theory NLCS is described as the spontaneous emission of one photon by the 
electron interacting with an intense external electromagnetic field. The electromagnetic field 
is described classically, an approximation which is justified for high intensity; the most used 
model until recently was the monochromatic plane wave. The theoretical studies published in 
the last three years adopt a more realistic description of short pulses, the pulsed plane wave 
model, for NLCS [9jfl3] as well for other processes [T^l H3] . The finite transverse extension 
of a real laser beam is neglected in quantum calculations, where an adequate treatment of 
this aspect was not developed up to now. On the contrary, in calculations based on CED, 
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where it is possible to include any shape of the electromagnetic field, beam size effects have 
been already considered [TBI E] • 

From a predicted probability distribution for the simultaneous detection of the emitted 
photon and the scattered electron, quantum theory extracts separate information on the 
emitted radiation, to which the majority of the results in the literature refer, or on the 
electron distribution, less studied up to now. It was remarked [10] that in experiments 
the radiation emitted in NLCS was not investigated, contrary to the Thomson case where 
the angular distribution of first several harmonics has been recorded [IB]. Electrons have 
been detected in the E-144 experiment at SLAC performed 15 years ago: the collision of 
a 46.6 GeV electron beam with terawatt pulses from a Nd:glass laser at 1054 and 527 nm 
wavelengths. Comparison with the theoretical energy spectra presented in Fig. 4 of [1] 
reveals the absorption of up to four laser photons. More than that, a suplementary evidence 
was given by the detection of positrons [5] which come from a succession of two elementary 
processes: NLCS, and (e~,e + ) pair creation by the energetic photon emitted in the first 
process. 

In the monochromatic case, the electron energy distribution presents thresholds [H [5] that 
can be connected with the absorption of fixed number of laser photons; these thresholds will 
be discussed in Sect. IIHI 

In connection with the perspective of very intense sources of radiation [19] , the interest of 
studying NLCS along with other processes in the very intense regime (intensities above 10 22 
W/cm 2 ) has increased and it is reflected in the most recent works [20] • Theoretical aspects 
as the use of a wave packet for the description of the initial electron [21] or the quantum 
description of the external electromagnetic field are reanalyzed [22J. 

Our paper is a theoretical study of the electron distribution in NLCS occurring in the 
interaction of the electron with a laser pulse. We have recently published a few results 
from a preliminary investigation of this distribution [23] • The present paper is organized as 
follows. Section [IT] displays the general expression of the multiple differential distribution 
describing both photons and electrons, from which analytic expressions for the electron 
distribution in the monochromatic or pulse case result. The monochromatic case is discussed 



in Sect. Ill Our purpose is not a numerical calculation of the transition rates, but an 
analysis of the position of the singularities they present in this ideal case. In the study of 
the electron distributions in the monochromatic case, one has to distinguish between the 
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"bare momentum" p 2 (the momentum of the asymptotically free electron) and the "dressed 
momentum" q 2 of the scattered electron; these momenta are in biunivocal correspondence. 
We have found that the analytic part of this exploration is more easily done in terms of the 
dressed momenta. The electron distribution present 5-type singularities at particular values 
of the final dressed momentum q 2 ^ , indexed by the positive integer N which is interpreted 
as the number of laser photons absorbed by the electron. We follow these singularities in 
the space of the polar coordinates ( | q 2 1 , 9 q2 , (p q2 ) of the vector q 2 . The manifold of the 
singularities for a fixed N is a surface which may be closed or not. We find a subset of 
points on this surface at which the distribution has a particular type of singularity that 
influences the electron angular distribution. Then, we translate the results in terms of the 
bare momentum p 2 of the electron and use them in the numerical illustration that concludes 



Sect. Ill We find that, for not too high laser intensity, this type of representation is useful 



in understanding the distributions obtained with finite pulses. Section [TV] is devoted to the 
equations valid in the pulse case. The numerical results presented in Sect. [Vjfor two type 
of pulses reveal conditions in which the analysis done in the monochromatic case is useful. 
In Appendices [A] and [B] we give details about the justification of some of the results 



presented in Sect. Ill concerning the monochromatic case. 



II. THE THEORETICAL FRAMEWORK. GENERAL EXPRESSION FOR THE 
TRANSITION PROBABILITY 

The general theoretical framework is the same as in [9], so not too many details will be 
given here. The same formalism is described in other recent publications [TUHT2] . 

The system investigated consists in an electron (charge e < 0, mass m ) interacting 
with the quantized electromagnetic field describing the emitted photon and with a classical 
electromagnetic plane wave with the unit vector ni in the direction of propagation and 
described by a vector potential 

A(0), = t_ ni .r/c=^, (1) 

where c is the velocity of the light and n\ = (1, ni) is the notation for a four vector 
associated to the laser propagation direction. We take the z axis along rix = e z and we 
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work with the vector potential 

A(0) = A f(4>) [ e x cos(C/2) sin(wi <t>) + e y sin(C/2) 003(0;! 0) ] . (2) 

This way the unit vectors e x and e y of the x and 7/ axes are chosen along the axes of the 
ellipse associated with the state of polarization of the external field. The wave polarization is 
described by the parameter ( ( £ = and tt correspond to linear polarization, £ = ±7r/2 to 
circular polarization); in our numerical calculations only the case of circular polarization will 
be considered. In the monochromatic case, f = 1 and U\ is the laser frequency; in the case 
of a laser pulse, the function / is the pulse envelope, supposed to be significantly different 
from zero only in a finite interval (4>- m , 0f), and uj\ is the central frequency. The maximum 
amplitude reached by the electric field of the pulse is ui\ Aq and the electromagnetic field 
intensity is characterized by the dimensionless parameter 

V = • 3 

mc 

The formalism we use starts with the general definition of the transition amplitude be- 
tween two states of the system (electron + quantized electromagnetic field + classical elec- 
tromagnetic plane wave), 

Mi_ a = lim lim (^ 2 (t 2 )\U(t 2 ,t 1 )\^ 1 (t 1 )) , (4) 

<2— >00 tl— >~ OO 

with U the evolution operator of the system. The initial and final states are products 
of free electron states of momenta pi and p 2 with, respectively, the vacuum state of the 
electromagnetic field and the one photon state of momentum k 2 and polarization S2. In 
contrast to [9j, were a spinor with well determined momentum, normalized on an arbitrary 
volume V , was used for the initial state of the electron, now we describe formally the 
electron with a momentum p! by a "wave-packet", 

*i(r, t) = J Hp f^ly/! e(p) dp , I *(p) | 2 = 6(p - Pl ) , (5) 

where £(p) is a solution of Dirac equation (p — mc) £(p) = and it is normalized to 1 
= 1), and E = ^m 2 c A + c 2 p 2 . This procedure was used recently in |2T] for a spinless 
particle. 

The interaction of the electron with the quantized electromagnetic field, responsible for 
the photon emission, is treated in the first order of perturbation theory. The action of the 



free evolution operator Uq (describing only the electron in the classical electromagnetic 
plane wave) on the free states leads to the Volkov states for which we use the explicit 
expression in Eq. (B4) of [PJ, with V replaced by (2nh) 3 . In the following formulas we 
use the four-momenta pi, p 2 of the initial, respectively final electron [pj = (Ej/c, pj) , 



Ej = cym 2 c 2 + p|, j = 1,2] and k 2 of the photon. 

The expression of the transition probability for the emission of a photon with the wave- 
vector k 2 G dk 2 and a scattered electron with momentum p 2 G dp 2 , averaged over the 
initial spin of the electron and summed over the final spin, is : 



d 4 n unpo i = n 4 (/c 2 ,p 2 ) 5(pi± - P2± - ^k 2± ) 6[m ■ (pi - p 2 - hk 2 )\ dk 2 dp 2 



(6) 



where the subscript _L is used to indicate the components orthogonal on the laser propaga- 
tion direction ni and the function II4 has the expression 



n 4 (/c 2 ,P2) 



el m 2 c 5 



1 



4 ir 2 E 2 (ni • pi) hw 2 



B? -1 + 



(pi • ni){p 2 ■ k 2 ) + (pi ■ fc 2 )(p 2 ■ ni) 



(M 2 ■ m) 



2(n x -pi)(7n -p 2 ) 



(ni ■ pi)(ni ■ p 2 ) 



mc) 2 k 2 ■ ni 
Api Ap 2 



+ 

, (7) 



(mc) hk 2 ■ ni 1 n i ' P2 

( eo = e/V^eo )• The external field dependence is contained in three one-dimensional 
integrals B, A x and A y defined as 



#(2,1) = 

A(2,l) = 
where the function G(pi,p 2 , fc 2 ; 0) is 



exp 



G(p!,p 2 ,k 2 ; 



eA(0) 



mc 



exp 



G(pi,p 2 ,k 2 ; 



(8) 
(9) 



G{p 1 ,p 2 ,k 



2: 



c 2 n i ' (Pi ~ P2 - ^2) + F(pi; 



F(p 2 



F(p;<f>) = —^- /rfx[e 2 A 2 ( X )-2eA( X )-p| 
2 ni • p J 



(10) 



(11) 



</><> 



with n\ 



-ni). In the case of a pulse, where the vector potential is different from for 



>in, W 



one has 



In the monochromatic case the indefinite integral can be 



used, as the change of the value given to 0o leads only to the modification of a phase factor 
in the Volkov solution. 
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Finally, we remind here a classicality criterion presented several times in the literature 
(see, for example, [T01 H3]) : the scattering of the radiation can be treated in the framework 
of CED, if the ratio ( s e g ftjjj\ 71 )/( m c 2 (1 + rf) ) , with 71 the Lorentz factor of the initial 
electron and s e ff the maximum number of laser photons absorbed, is small compared to 1. 
For 77 < 1, s e fr is of the order of unity; for rj 1 , s e s increases rapidly, proportional to 
rj 3 , and the ratio becomes 

y = ^. (12) 

mc z 

If y becomes of the order of unity or larger, then the quantum behaviour sets in, and, as 
discussed before, one can expect to obtain a final electron distribution different from the 
initial one. 



III. THE MONOCHROMATIC CASE 

In the monochromatic known for long time [T], the integrals A x ,A y and B 

have analytic expressions as series of generalized Bessel functions. In these series each term 

contains an one-dimensional 5-function, as illustrated here by the integral A x , 

00 

•4,(2,1)= Ai N ^[n 1 -(q l + Nhk 1 -q 2 -hk 2 )]. (13) 

N=-oo 

The four-momentum k\ , 

k 1 = (k° 1 ,k 1 ) = —n 1 , (14) 
c 

is interpreted as the momentum of a photon associated to the electromagnetic monochro- 



matic plane wave. In Eq. (13) appears the dressed four-momentum q, a quantity met 
also in the description of the electron motion in classical theory, connected with the bare 
four-momentum p by 

m\Jp TT e 2 Al . . 

q = pA m, U P = - , n l -p = n 1 -q. (15) 

rii ■ p Am 

Between the two 4-momenta q and p the correspondence is biunivocal, p can be expressed 
as a function of q as 

p = q n\. (lb) 

m ■ q 

For a given four momentum p, the first component is po = E/c and the first component of 
q is go = W/c, where 



7 e 2 J[2 

E = cJ p 2 + m 2 c 2 , W = cJ q 2 + m 2 c 2 , m* = m\l + - — ^ , (17) 

2 m z & 



with m* the dressed-mass, named also the shifted mass. 

The use in the fully differential probability ^ of the integrals -4(2, 1) and 13(2, 1) , 



as series of S- functions similar to (13), gives an expression for <i 4 n unpo i from which, after 
standard manipulations, one extracts the transition rate, denoted d 4 T. It has the structure 
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d 4 T(p 2 , q 2 ) = 7f } (q 2 , k 2 ) 5( qi + Nhh -q 2 - kk 2 ) dq 2 dk 2 , (18) 

N=l 

i.e., d 4 T is a series of four-dimensional 5-functions with coefficients depending on the vari- 
ables q 2 and k 2 . A term with fixed iV in the previous expression is the contribution to 
the differential rate of the process in which N laser photons have been absorbed. 

In the monochromatic case it is customary [21] to present these distributions as functions 
of q 2 , but it is also possible to present them taking as variable p 2 , using the relation 

/ mcUp \ , 
dq 2 = 1 + f— dp 2 , 19 



with the ponderomotive potential Up defined in (15). 

In the following we shall suppose that in the "partial rates" 74 the connection be- 
tween the momenta imposed by the 5-function was observed. In fact the product of four 
one-dimensional (^-functions leaves arbitrary only two of the six components of the three- 
dimensional final momenta q 2 and k 2 . Our purpose is not the evaluation of the partial 
rates, but the analysis of the implication of the conservation rules for each term with fixed N, 
in the case of the electron distribution. As we shall see in Sect. |Vj in appropriate conditions, 
connections are possible between the results of this analysis and the electron distributions 
in the pulsed wave case. 

In order to get the differential rates describing the electron energy and angular distribu- 
tions we have to integrate over the emitted photon momentum. This is a direct operation 
performed by three of the four 5-functions in each term and it imposes the following values 
to the emitted photon momentum 

hk 2 = qi - q 2 + TVfrki . (20) 

The expression allowed for the frequency is ck\= c \ k 2 | . After the integration on k 2 only 
an one-dimensional 5-function is left in each term of the series which represents the double 
differential rate d 2 T e describing the scattered electron, in terms of the dressed momentum 
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Q2, 



d 2 r e = ^r^(q 2 )5 



N=l 



1 



m*c 



■(q° + hk°-q°-Nhk[) 



dq 2 



(21) 



with r( JV )(q 2 ) = T f ) (q 2 ,k 2 )/(m*cft 3 ). 

In the following we emphasize some particularities of the electron distributions that come 



out from an analysis of the argument of the 5-function in (21). 



A. Simultaneous detection of electron energy and direction 



In the monochromatic case the final electron distribution is written as a series of 5- 



functions, as displayed by Eq. (21). We shall study the position of the singularities in the 



variables W2 [ or, equivalently, | q 2 1 related to W 2 by (17) ], and q 2 , the unit vector along 
the direction of the final dressed momentum q 2 of polar angles 8 q2 and (j) q2 . It is convenient 
to introduce a new four-vector 



Qn = <?i + NHkt 
for which we have 



,£ 



N 



Q N ), Q N = qi + N hki 



Q% = (m*c) 2 + 2Nhk 1 - q x 



£ N = W 1 +N hw x 



(22) 



(23) 



We emphasize that for given laser intensity and fixed N, the four- vector Qn is well deter- 
mined only by the values of qi and k± . In the reference frame described at the beginning 
of Sect, [nj we denote by 9^ and <pN the polar angles of the vector Q N , and by «jv 6 [0, ir] 



the angle between the momentum q 2 and the vector Qjy • With (22) the argument of the 



5-function in a term with fixed N in (21 ) is 

1 



The equation 



iq o + hk o_ Q N) 



;(W 2 + c\ Q N -q 2 I -£. 



N 



(24) 



F N {q 2 ) = 0. (25) 
determines the position of the singularities in the space of the variables W 2 (or | q 2 | ) and 



q 2 . The condition (25), for N > 1 defines a family of surfaces in the space (|q2|, 0. 



the differential rate (21) has a 5-type singularity along these surfaces and is zero otherwise. 



Before going in more detail, we draw attention to an approximate symmetry property of 
-F/v , valid in the case 77 ~ 1, when the number of terms which gives practically non- negligible 



contribution to the electron distribution is limited to a value N max of the order of unity. In 
this case it makes sense to analyze the condition (25) only for N < N max . If, in addition, 
71 3> 1 and the direction of the bare momentum pi is not too close to n l5 then, the angle 
«at between q2 and Qtv is approximately equal to the angle between q2 and pi. As a 



consequence, the solutions of Eq. (25) have, with a very good approximation, a rotational 
symmetry with respect to the direction of the incident electron direction. 

We remark also that because in the regime (71 ^> i] ~ 1) the dressed momentum q 2 is 
very close to p 2 , the two distributions, one expressed in terms of the variables of q 2 , the 
other in terms of p 2 , are almost identical. 

In order to obtain the energy distribution or the angular distribution of the scattered 



electrons, one needs to integrate the double differential distribution (21 ) over the parameters 
that are not observed, by writing the (^-function in a way convenient for each distribution. 
For the angular distribution we need the relation: 

8{X - X so1 ^ 

sol 



5(F N ( q2 )) = 



dF N 
dX 



m*c 



sol 



Q°2 


Qn - Q2 




\Qn 


y/cOS 2 OtN - 


- Cn 



6{X-X 



sol\ 



(26) 



which is based on the solutions (A9) and (A10) of Eq. (25) for the unknown X = |q 2 |/(m*c); 



the quantity Cn is defined in (A2). By the generic summation index 'sol', we understand 



the (one or two) solutions acceptable at fixed N (see details in Appendix A). When this 



expression is replaced in (21), it displays the position of the singularities in X at fixed 
direction of q 2 . 



For ^i^- = the expression in (26) for <5(-Fjv(q 2 )) is not valid. In the following, we 



dx 



work with cosajy < VCn and after obtaining the angular distributions, we take the limit 
cos«7v = VCn- As shown further, the singularity present in the double differential distri- 



bution (21) influences the angular distribution of the electrons, obtained after integration 



on the scattered electron energy. 



To prepare (21) for the calculation of the energy distribution, we have to find the polar 



angles of q 2 that are solutions of (25) at fixed | q 2 | . The equation (A4) gives us, for any 
X in the interval (A5), the unique solution for cosaAr, denoted by Y(X) (see Appendix 
A). From it we derive the possible values for the polar angles of q 2 by solving the equation 

cos 9 at cos 6 q2 + sin N sin 6 q2 cos(0 g2 — 4>n) = Y(X) , (27) 

considering as the unknown variables one of the two angle 6 q2 or (fi q2 , with fixed Y(X). 
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Simple particular cases are collinear and head-on collisions [sin^jv = 
when Eq. (27) is an equation for 9 q2 only, with the solution cos 9 q2 



, cos 9 jy = ±1 ], 
= a N Y(X) with 



<j n = sgn(cos 6 N ). In this case we write the ^-function in (21) as 

IQjv - O2I 



S(F N (q 2 )) = m*c 



5(cos6 q2 — ai<[Y(X)) , sin^jy = . 



(28) 



lOal IQivl 

The cases 9 N ^ 0,tt are more complicated since both angles 9 q2 and (f> q2 appear in Eq. 
(27). One possibility is to solve Eq. (27) for the unknown <fi q2 , keeping as parameter 9 q2 . 



As shown in Appendix B, Eq. (27) has two solutions, 

Y(X) — cos 9n cos 9 q2 



. sol 



± $0) "^0 = arccos 



sin 9 N sin 9, 



'12 



if 9 q2 obeys the condition 



cos 9 q2 G [cos(9n + «jv), cos(9n — on)] ■ 



(29) 



(30) 



This condition defines an angular range that we denote by Xg. We emphasize that in 
the present context the value taken by the angle a* depends on |q2|, being expressed as 
otN = arccos(F(X)). Finally, the procedure leads to the expression of ^-function in (21) 

. sol\ 



KFnM) 



m*c 



sol 



\Qn 


-02 


1 


I02I 


\Qn\ 




sin 9n sin 9 q2 sin $ 



sin^^O. (31) 



With this expression of 8{Fn) the distribution (21) displays the position of the singularities 



in (p q2 at fixed 9 q2 and | q 2 



If we choose to solve the equation (27) for the unknown 9 q2 with cf) q2 as a parameter, the 
solutions are more complicated. With the notation u q2 = cos9 q2 , one finds (for details, see 
Appendix B) two possible solutions: 



u §° = — ( Y ( X ) cos ^ ± sin 0*1 cos(0 92 - <f> N )\ ^s N - Y 2 (X)) . 

S I\T V / 



S N 

sn = cos 2 9 n + sin 2 9n cos 2 



''12 



>N 



(32) 
(33) 



Depending on the initial conditions and on the value of (p q2 , one or both solutions are 
acceptable, namely: 

i) for Y 2 (X) < cos 2 9n, only one solution is acceptable for any value of <p q2 : u& , if 



cos9n cos 



''12 



< 0, and u q2 , if cos9n cos(4> q: 



> 



ii) for cos 9n < Y (X) < 1, both solutions acceptable, but the domain of cf) q2 is 



e 7r- 



reduced to <fi q2 G [4>n — 0o, (pN + 4>o] if Y(X) > and to 
if Y(X) < , where O = arccos \/(Y 2 (X) - cos 2 a N )/ sin 2 a* G (0, 7r/2) 



i> O ,7T + 0jV + 0O. 
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The 5-f unction in (21) is written now as 



5(F N (\ q 2 \;u q2 ,(f) q2 )) = m*c > . ; ; — : 

z H4jv| |q 2 | 



52 M <?2 



(sol) 



It gives the position of the singularities in 6 Q2 at fixed (p q2 and | q 2 |. 



(34) 



B. Angular distribution of electrons 



We get the angular distribution of electrons using the expression (26) in the distribution 

(ra*c) X, with the result 



(21) and integrating over |q 2 

dr. ^ 

N>1 sol 



dtt 



Q2 





2 W 2 


IQjv - qa| 


c 







rW(q 2 ) 



Vcos 2 «at — C. 



iV 



(V) 
e,ang' 



(35) 



|q 2 |=m„cX so1 

The substitution rule indicated above means that the modulus of q 2 must be replaced 



everywhere by (m*c)X so1 , with X so1 given by Eq. (A9) or (A10). 

Based on the results in Sect. HI A[ Appendix [A] and some more details given in Appendix 
[B] we mention here the main features of the angular distribution. We describe the situation 
of a term T^a ng , with a fixed value of N. If, for that N, we are in the case I, when < 0, 



[Cn defined in (A2)] there is one solution A so1 [Eq. (A9)] for any cosa/v G [0,7r], i.e. 
for any direction of the scattered electron, and the term ri N J n g in the sum (35) is finite. 
If, for the considered N , we are in the case II, when < Cn < 1, then, there are two 
solutions A so1 [Eq. (A10)] for any direction o>n obeying the condition cosaTv < s/Cjy. For 
cosojtv = VC-n, the two solutions X± in (A10) coalesce and the corresponding T^ang has a 



singularity. The condition cosaTv < VCn, determining the possible scattering angles for a 
given N in the case II, can be expressed in terms of polar angles of the electron in the form 
{9q 2 ,(fiq 2 } £ V{Cn)- The explicit expression of the domain V{Cn) is deduced in Appendix 
B. 



C. Energy distribution of electrons 



The energy distribution is obtained by integrating the fully differential distribution (21) 
over the electron directions determined by the angles 6 q2 and 4> q2 . 
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For collinear and head-on collisions, using (28), the integral over 9 q2 is performed directly 
and the energy distribution becomes 

2vr 



dWo 



c 



E 

7V>1 



W 2 \Q 



N 



02 1 



'92 



IQ 



q2 



A' 



<?2 — s °l 



iV>l 



where so i = arccos(F(|q2|/(m*c))) , if = 1, and 9 so \ - 
un = —1 , with Y(\q 2 \/(m*c)) calculated according to (A4). 



J2 r i N w> sm6 N = 0, (36) 
7r — arccos (F ( | q2 1 /(m*c))) if 



For other initial configurations, it is convenient to use in (21) the expression (31) of the 
^-function and, as a consequence, in the calculation of the energy distribution the integral 



over 



J H2 



is performed directly After that, for the integral on 9 q2 that has to be done 



numerically, the domain of integration reduces to the interval Tq defined by the condition 

W 2 \Q N ~ qal 



(30). The final result reads: 
dT e in* x — < \ - 

N>1 sol 



de, 



<I2 



dW 2 

with so1 and $ given by (29) 



sin^Arsin^r 



IQ 



Q2j 



N\ 



N>1 



(TV) 
e,W 



(37) 



D. An example 

We illustrate the previous analysis by an example. We have seen that the <5-f unction in 



the multiple differential distribution (18) imposes the restriction (25) on the vector q 2 and 



we have described the position of the singularities in terms of ( W 2 , 9 q2 , 4> q2 ). As mentioned 



in Sect. Ill A the analysis can be converted in terms of the bare momentum p 2 . In this 
case we think of surfaces in the space ( E 2 , 9 P2 , (f) P2 ) on which the singularities are localized. 
In the example that follows we shall present graphs with the curves giving the position of 
the singularities in the plane {E 2 , 9 P2 ) at fixed <f) P2 — 4>pi- 

We choose the case of an electron of energy E\ = 46.6 GeV scattered by a circularly 
polarized monochromatic wave with the frequency oj\ = 0.043 a.u. (1.17 eV) and the field 
intensity / = 4.4 x 10 17 W/cm 2 (rj = 0.6); these conditions are close to those in the SLAC 
experiment, in which the detection of NLCS was achieved. We consider two cases for the 
initial direction of the electron: (a) 9 Pl = 0.9 ir, close to the value used at SLAC, and (b) 
9 Pl = 0.5 7r (orthogonal geometry), with Pl = in both cases. For these initial conditions 
we have n x ■ qi > m^c and we are in the case II (defined in III. A) for values of up to 
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pa 10 6 ; this value is much larger than the maximum value of iV contributing to the electron 
distribution at the intensity considered, which is of the order of ten. 

We describe the position of the 5-type singularities in the double differential distribution 



(21) based on Eq. (34), giving some details valid in our particular case, <f> Pl = (p P2 = 0, 
using as variables the bare energy E P2 and polar angle 9 P2 . In the present discussion, 
preceding Fig. 1, we have in mind only low values of N (of the order of ten), for which, as 
we have mentioned before, we are in the case I. In the particular case (f> P2 = <p pi = we have 



chosen, for any X in the interval (A5) (i.e. for any energy W 2 in (A8)) the two acceptable 



solutions of Eq. (27), given by Eq. (32), reduce to 9 q2 = 8 N ± arccos F(X) which coalesce 
for Y(X) = 1, i.e. at the ends W 2 = W a and W 2 = Wb of the interval. The maximum 
domain of variation for the angles is given by the condition 

q2 E [ 9 q ,A = 9 N - arccos F m i n , 9 q<B = 9 N + arccos Y min } , (38) 



with given by (A6). This domain can be transcribed in terms of bare energy and 



scattering angle, using the relation 



, |q 2 |cos^A,B -m 2 cV/( 4 ™i ^2) . ,, )(V . 
9 q A b — ^ 9 p A b = arccos (39) 

\J{W 2 — m 2 c^if /(Arii ■ q 2 )) 2 /c 2 — m 2 c 2 



for the angles and Eq. (17) for the energy. 

In the regime discussed here ( 7 ^> rj ~ 1 ), qi and q 2 are very close to the corresponding 
bare momenta pi and respectively p 2 , so the results in the plane (E P2 ,9 P2 ) are practically 
identical at the graphical level to those in (W 2 , 9 q2 ); in particular, we have W a j, ~ E a ^ 
and 9 p a,b ~ 9 q a,b- Another particularity is that the upper limits of the intervals Wb(N) ~ 



Eb(N), defined in (A8), are almost independent of and approximately equal to the initial 
electron energy E\\ the lower limits, however, are significantly dependent of N . Then, the 
energy of the final electron in the process in which iV photons are absorbed takes values in 
an interval E 2 € (E a (N), Ei) with E a (N) < E a (N — 1). The successive values E a (N) are 
named thresholds of the energy spectrum. 

In Figures [l] (a) and (b), for the two values of 9 pi we have chosen for the direction of the 
incident electron, we display in the plane (E P2 ,9 P2 ) the curves that represent the solutions 



of Eq. (25) for iV taking values from 1 to 11. The coordinate along the y axis is 59, 
defined as 69 = 9 P2 — 9 Pl ; the good symmetry of the two figures with respect to the value 
69 = is a consequence of the rotational symmetry around the direction of pi mentioned 
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FIG. 1: Position of the first eleven lines in the plane (E 2 ,9 2 ) f° r E\ = 46.6 GeV, rj = 0.6, circular 
polarization and two directions of the incident electron: (a) 6 Pl = 0.9ir and (b) 8 Pl = 0.5tt; (c): 
detailed view of (a) for N = 1 . 



in Sect. HI A The results show that 89 takes very small values for all N, i.e. the final 
electron direction is very close to the initial one for all the cases represented (N < 11). On 
the other hand, the energies E 2 are spread till relatively small values. The interpretation 
of these results is that in case (a) the electron can lose up to 45% of its initial energy in 
the process in which only one photon is absorbed (N = 1), and up to 85% for N — 11. 
In case (b), when the initial electron incident orthogonal on the laser pulse direction, the 
angular distribution is more compressed towards small angles and it is also compressed in 
the region of larger final energies. As, according to the conservation laws, the energy lost 
by the electron is converted in the energy of the emitted photon, this means that this case 
is less efficient for energy conversion. 

In Fig. [I] (c) is presented only the curve with N = 1 from the ). The limits 

E a (l ) and -£^(1) of the domain in energy that gives contribution to the spectrum, obtained 



from Eq. (A8) using Eq. (17), are marked on the graph; for any E 2 G (E a (l), -E'b(l)) there 



are two angles 9 q2 for which the equation (25) is verified, they become a double root for 
E 2 = E a (l) and E 2 = E b {l). As discussed after Eq. (39) the value E a {l) corresponds to 



the threshold of one photon absorption in the energy spectrum. The domain of angles 6 P2 
that contribute to the angular distribution for N — 1 is the interval (9 Pl + 89i, 9 Pl + 89 2 ); 
within this interval there are two solutions E 2 of Eq. (25), which coalesce for 9 P2 = 9 Pl +59± 
or 9 r , 



9,. 



59*. 



'P2 ~ "Pi 1 wJ 2- 

Another aspect worth to be discussed is which would be the dependence of the type of 
curves represented in Fig. [I] (a) and (b) on the laser intensity. One feature to be considered 
is the increase of the maximum number of photons that can be absorbed in the process with 
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the laser intensity. The other feature is that, at fixed N and q±, when the laser intensity 
increases the curves tend to become closer to each other, i.e. the domain in which the 
energy of the final electron can take values shrinks. This process can be understood using 
the concept of "dressed mass": when the laser intensity increases the electron becomes 
heavier, and consequently its recoil at fixed number of photons absorbed decreases. 



IV. THE PLANE WAVE PULSE 



We consider now the more realistic model of a pulsed plane wave, going back to Eqs. 
(|6| and ([7]). In the pulse case, the integral B is expressed in terms of convergent integrals, 
using Eq.(30) of [9] (see also, [H] for an alternative approach). 

Now, as only three (^-functions appear in ([6]), only three conditions are imposed to the 
six variables p 2 and q 2 , namely 

Pi± - P2± - frk 2 _L = , m • (pi - p 2 - hk 2 ) = . (40) 

In order to get the one-particle (electron of photon) distribution, the differential distribution 
(|6| is integrated over the momentum of the that is not detected using the conservation rules 



(40). The integration over the orthogonal components of the momenta is performed directly 



using the 5 function, so from (40) we get the replacement rules to be used in ([7]), 

P2± -»> P2± = Pi± - Kk 2 ± (41) 

for the integration over the orthogonal component of the photon momentum, and respec- 
tively, 

hk 2 ± -)> hk 2 ± = Pi_l - P2_l (42) 

for the integration over the orthogonal component of the electron momentum. The integra- 
tion over the third component requires some further calculation, due to the fact that the 



second relation in (40) contains a combination of energy and momenta. We present the 



results in both cases. Using the adequate relations from the following ones, 

S[n ■ (pi -p 2 - hk 2 )\ = — rr\ 5 G°2z -p~2z) = r 5 (k 2z - k 2z ) , (43) 

cm ■ (pi - hk 2 ) cm ■ (pi - p 2 ) 

with 

~ _ (mc) 2 + (pii ~ hk 2± ) 2 _ m ■ (pi - hk 2 ) ^~ (pi± - p 2 ±) 2 _ n x ■ (pi - p 2 ) 
2n l -{p 1 -hk 2 ) 2 " 2n 1 -{pi-p 2 ) 2 



(44) 



16 



and 



E2 = imcf + ( Pl± 
c 



hk 



2± 



Tlx ■ [pi - Hk 2 



frjjJo 



(Pll - P2±) 2 , ni ■ (pi -p 2 ) 



2ni ■ (pi - ftfo) 2 c 2ni ■ (pi - p 2 ) 

one obtains the two one-particle (photon or electron) distributions, 

d 2 I\ = 



E ui 2 

— ^Il 4 (k 2 ,p 2 )duj 2 dtt k2 , 



d 2 T P 



cni ■ p 2 c 
E 2 uj 2 I P2 



cni • ^2 ^ 2c2 
The attribute unpolarized was omitted. 



n 4 (&2,P2) rf^ P2 • 



(45) 

(46) 
(47) 



NB. The quantity denoted here by k 2 is different from that defined in (20) and used in 



Sect. |III[ as it comes out from a different conservation rule. 

The structure of the previous two distributions, using each in a specific way the same 
function II4 implies the possibility of connections between the two distributions, as it will 
be mentioned at the end of Sect. IV Bl 



V. NUMERICAL RESULTS 



We consider two type of pulses: i) a pulse with a finite duration, of almost rectangular 
shape, with the envelope /(</>) in Eq. ^ constant on a region of length equal to a multiple 
N c of periods of the carrier, and two very short smooth wings; we shall name this pulse 
"quasimonochromatic" , ii) a pulse without a constant region, consisting in two wings of 
variable length. We shall see that the first type of pulse, if N c is large enough, leads 
to results similar to those predicted by the monochromatic approximation, which explain 
the adopted terminology. The contribution to the scattering probability of the wings of 
the pulse is very small compared to the contribution of the flat central region, still these 
smooth wings are required in order to ensure the continuity of the vector potential and of 
its derivative. The results obtained with the second pulse are considerably different from 
the monochromatic ones. 

In Fig. [2] is represented the vector potential A as a function on <f>/T for the two pulses 
mentioned before. The rectangular portion of the pulse in (a) has the length of N c =8 
periods of the carrier. In the case (b) we have chosen a cos 2 envelope, the total length of 
the pulse corresponds to a number N t of periods of the carrier equal to 10. In both cases 
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FIG. 2: (Color online) The components A x and A y of the vector potential for a rectangular pulse 
with N c = 8 cycles (a) and for a cos 2 pulse with N t = 10 cycles. 

the parameter rj is 0.6 ( the value used in the SLAC experiment). In all numerical examples 
presented here we choose the laser central frequency u>i = 0.043 a.u. (1.17 eV) and circular 
polarization. 



A. Effect of the pulse shape 

We have calculated the electron double differential probability distribution d 2 T e /dE 2 dQ P2 
for the conditions of the SLAC experiment (u\ = 0.043 a.u., rj = 0.6, E\ = 46.6 GeV) and 
for the two pulses represented in Fig. [2} 

In Fig. [3] (a) and (b), we present results in a logarithmic color scale in the same coor- 
dinates as those in Fig. [I] (E 2 , 56 = 6 P2 — 6 Pl ). The case considered is 4> P2 = 4> P1 = 0, as 
111 Fig. [1} Due to the very good symmetry of the results with respect to 56 = 0, remarked 



also in the discussion of the monochromatic case (Sect HID), only the values 56 < are 
presented. In both figures one can see a series of maxima located on curves with the same 
shape as those presented in Fig. [T] in the monochromatic case. Notice that in Fig. 1 both 
56 < and 56 > are represented. For the intensity we consider, with a relatively low value, 
only the first 11 maxima are visible. The difference between the two cases is that, while for 
the rectangular pulse the maxima are very sharp, and have a fine substructure, for the cos 2 
pulse the main maxima and their subpeaks become wider and smooth. 

In Fig. [3] (c) are represented the energy distributions dT e / dE 2 , for the rectangular pulse 
(full black line) and for the cos 2 pulse (dashed red line). The two distributions are similar, 
having a "ladder-like" structure, with successive shoulders which can be understood based 
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FIG. 3: (Color online) (a): Electron probability distribution ggr^f — for the conditions of SLAC 
experiment and for a rectangular pulse; (b): the same as in (a), but for a cos 2 pulse; (c): energy 
probability distribution for the two pulses; (d): angular probability distribution for the two pulses. 
Black (full) lines refer to the case of rectangular pulse, and red (dashed) lines are for the cos 2 pulse. 

on the monochromatic limit: their positions coincide with the thresholds E a (N) « W a (N) 



in Eq. (A8); the upper limits Ef,(N) are almost independent of A" and approximately equal 
to the initial electron energy E 1 . The first interval at the right in Fig. [3] (c) covers the 
region E 2 E (E a (l) ps 27.6 GeV, E\) and can be interpreted as the sum of contributions of 
the processes in which any number of photons N > 1 can be absorbed. The next step, the 
region E 2 G (E a (2) ~ 18. GeV, E a (l) 27.6 GeV), is the contribution of the processes with 
A^ > 2, as N = 1 does not contribute anymore and so on. The fact that the values of the 
successive steps decrease very fast (note the logarithmic scale) is due to the relatively small 
value of 7], still close to the perturbative regime. The figure [3] (c) is similar to Fig. 4 in [3], 
calculated there in the monochromatic approximation. 

The angular distribution dr e /dfl P2 for <p P2 = <p Pl = and variable 9 P2 for the same two 



pulses as before is presented in Fig. |3| (d). Here one can see again the same "ladder-like" 
structure, but, unlike in the case of the energy distribution, there is a sharp maximum at 
the left end of each step. These maxima are the corespondent of the singularities of the 
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FIG. 4: (Color online) Electron probability distribution 



d 2 r c 



dE2dQ p 



as a function of 72 and 56 for head- 



on collision, initial electron energy E\ = 10 mc 2 , and four values of the field intensity: r\ = 0.5 
(a), 77 = 1 (b), 77 = 2 (c), 77 = 5 (d). 



angular distribution existent in the monochromatic case [see Sect. Ill C ] and are localized 



at 59 p a(N) given by (39), as presented in the example III.D. As expected, they are much 
better defined for the rectangular pulse than for the cos 2 pulse. 



B. Effect of the field intensity 

We illustrate now the influence of the laser intensity on the double differential distribution 



(47) of the electron. 

In Fig. [4] we consider the case of a rectangular pulse with N c = QT and a head-on collision 
(Opi =7r) with the electron having the initial Lorentz factor 71 = 10 4 . As in the previous 
subsection, we represent the probability distribution in the plane (72, 56 = P2 — 9 Pl ). 

We remind that since we are in the case of head-on collisions and the laser is circularly 
polarized, the electron spectrum is symmetric with respect to rotations around the z axis, 
i.e. it does not depend on the angle 4> P2 . We have chosen four values of the parameter r\: 0.5, 
1, 2 and 4. For the first two values the spectra present a series of maxima localized along 
curves whose shape and distribution is that of the lines predicted in the monochromatic case 
(see Fig. [I]). However, when r] increases, the successive peaks become so close to each other 
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that they start to overlap, tending to form a smooth continuum, so we hardly distinguish 
them in Figj4] (b) and not at all in Figs. W(c) and (d). This behaviour is in agreement with 



the discussion at the end of Sect. Ill In the last two cases we remark another interesting 
feature: for 77 > 2 the distribution does nor cover uniformly the plane (72, SO) but only a 
small region, with a well defined shape. 

The behaviour found in Figs. [4] (c) and (d) can be understood in correlation with the 
photon distribution. In [25] it was shown that in the case of an ultrarelativistic electron 
and for large values of rj, the photons are emitted only in a well defined, very small domain 
of angles; although in the cited paper only the CED formalism is used, it can be shown 
that the conclusion concerning the photon distribution are valid also in the quantum case. 
For a rectangular pulse, as that considered here, and for head-on collisions, the emitted 
radiation has a continuous spectrum, extended from u 2 = and up to a maximum value 
Qm, an d it is emitted practically at a constant angle 8 , symmetrically around the z axis, 
@k 2 ^ {^0 ~ S,0q + 5), 5^1. In terms of photon momentum, this means that the func- 
tion 11^2, p 2 ) m ^ is non-negligible only for k 2 along the directions of the unit vectors 
eo = (sin #0 cos 0, sin 6q sin 0, cos 6q), with G (0,27r); when expressed in terms of electron 



momentum, according to the conservation rules (42), (45), this condition leads to the partic- 
ular shapes present in Fig. [1] (c) and (d). The correlation between the electron and photon 
distributions needs further investigation. 

C. Effect of the initial electron energy 



We present in Fig. |5jthe electron energy distribution dT/dE 2 for the case of a head-on 
collision, rj = 5, and for four values of the initial electron energy: ji = 10 (Ei = 5.11 
MeV) in full line, 71 = 10 2 (E 1 = 51.1 MeV) in dashed line, 71 = 10 3 {Ex = 511 MeV) in 
dotted line, 71 = 10 4 (Ei = 5.11 Gev) in dash-dotted line. The laser pulse is rectangular, 



with N c = 10 cycles. The values of the classicality parameter (12) in the four cases are, 
respectively, y = l.lx 10~ 4 , 1.1 x 10" 3 , 1.1 x 10" 2 , 1.1 x 10 -1 . 

The coordinate along the x axis is chosen as 1 — E 2 /Ei and the results are presented in 
a log- log scale. For the lowest value considered for the incident electron energy, the energy 
distribution decreases very fast with the ratio EijE\: practically the entire distribution is 
contained in the interval E 2 G (0.992£i, E\) which is an indication of the classicality of the 
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FIG. 5: (Color online) The energy distribution dT e /dE2 for r] = 5 and four values of the initial 
electron energy: E 1 = 5.1 MeV (full black line), E\ = 51.1 MeV (dotted red line), E 1 = 511 MeV 
(dashed green line), E\ = 5.1 GeV (dot-dashed blue line) . 

process. With the increase of the initial electron energy, the energy spectrum is extended 
towards lower values of E2/E1 , and its slope is much lower, as an indication of the onset 
of the quantum behaviour; for the largest value of the incident electron energy considered 
(71 = 10 4 , Ei = 5.11 GeV) the electron can lose up to 90% of its energy. 

The small peaks visible in the energy distribution for 71 = 10 3 and 71 = 10 4 correspond 



to the limit E a (l) « W a (l) [Eq. (A8)] of the energy range in which one photon absorption 
contributes. For the other two values of 71 these points are located at values of 1 — 
E2/E1 < 10~ 4 , not represented in our figure. The presence of this one photon peak can be 
explained using the analogy with the electron behaviour in the monochromatic case, based 
on conservation laws valid in that case: for E 2 very close to E\ the electron direction 
changes very little, and the photon is emitted at an angle extremely small with respect to 
the initial electron direction, i.e. 6k 2 ~ tt. It is known for a long time that for this 
geometry the terms with high N in the radiation spectrum are suppressed, i.e. only the first 
few terms contribute to the total rate, even if rj is relatively large, and successive thresholds 
are visible in the electron distributions. For smaller values of E2 the angular distribution 
widens, and many values of N contribute to the total rate, i.e. the typical behaviour for 
large rj sets in: the successive maxima become broader and overlap, giving rise to a smooth 
continuum. 
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VI. CONCLUSIONS 



The aim of our theoretical study of NLCS was a first description of several features of 
the scattered electron distributions in the case of a pulsed electromagnetic wave. We have 
identified two possible guides for a qualitative understanding of the electron distribution: 
the monochromatic limit and the emitted radiation pattern. Which one is useful, if any, 
depends on several parameters: pulse shape (duration and intensity) and initial electron 
momentum. While the role and condition of applicability of the first guide was identified 
in the present study, the connection with the emitted radiation spectrum requires further 
investigation. 



Appendix A: The solutions of equation (25) 



In order to present the properties of the solutions of the equation (|25|), we use the nota- 

Ulc 

Qjv| 



tions (22) and define a set of dimensionless quantities 

Q°n 



UN 



v N 



m*c 



and we use the notations 



cos oin = Y G [—1, 1] , C 



N 



X 



m\ 



> o 



(Al) 



it 



N 



v 2 N + l) 2 



4< 



(A2) 



It is useful to emphasize that while un, vn and Cn are determined by the initial conditions, 
X and Y are variables, connected with the final momentum q2 , which is subject to the 



condition (25). 



Analyzing the equation (25), we have found that we have to distinguish between two 
cases: 

Case I: Un — Ujy > 1, when Cn < 0, 
Case II: un — v n < 1, when < Cn < 1 • 

We have found that if the momentum of the incident electron fulfills the condition ni ■ 
<7i — m*c > 0, we are in the case I for Nhu\ > m*c w ^~ m * c2 anc l i n the case II for 



Wi—m*c 2 
ni-qi— m*c 



If n\ ■ qi — m*c < 0, we are in the case I for any value of N. 



The numerical examples presented in this paper (Sect. HID and IV) refer to the case 
of head-on or nearly head-on collision of a very energetic electron with a laser pulse of 
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moderate intensity (71 ^> i] ~ 1). For these initial conditions we have n\-q\ > m*c, and we 
are in the case II for all values of iV which gives non-negligible contribution to the electron 
distribution. 

When written explicitly, the function Fn [Eq. (p4|)] is a function of the two unknowns 



X and Y, defined in (Al) and (A2), respectively. This way Eq. (25) becomes 



F N (X, Y) = VX 2 + 1 - JX 



N 



2XYv N -u N = 0. 



(A3) 



It can be easily solved as an equation for Y, with the parameter X, leading to the expression 



Y(X) 



2u N y/X 2 + 1 - (u 2 N + 1 - 



N ) 



(A4) 



2Xv N 

The properties of the function Y(X) are different in the two cases mentioned before: 

Case I ( un > vn + 1): Y increases monotonously with X and the condition |Y(X)| < 1 
leads to a domain of acceptable values of X 



X e 



UN - V N 



2 2(u N - v N ) 

Case II ( Mat < t>Ar + 1 ): Y has a minimum 



u N + v N 



2{u N + v N ) 



(A5) 



reached for 



= VC N e (0, 1) 



^4u% -(u 2 N + l- v 2 N ) 2 
u 2 N -v 2 N + l 



(A6) 



(A7) 



The condition |Y(X)| < 1 leads to the same domain (A5) of values for X as in the 



case I. When expressed in terms of energy of the dressed electrons, the interval (A5) is 
W 2 e [W a , W b \ with 



W a {N) =£ N - 



Sn/c- I Q 



ni • gi, W b (N) = £_ 



N 



£n/c+ I Q 



ni • qi . 



N 



(A8) 



Going the other way around, i.e. solving the equation (25) for the unknown X as function 
of Y, we find: 

i) only one solution in the case I, namely 



X+{Y) 



v N Y(u 2 N -v 2 N + l) + 2u N v N VY 2 ~ C_ 



N 



2« - v%Y 2 ) 



(A9) 



for any Y G [—1, 1], 
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ii) two values for X in case II, 



Y / V x ^ y (tgr - < + 1) ± 2 Wjv yr 2 - c N 
±{ ) ~ 2K-<n (A10) 



for V G [K niT i, 1]. At Y = Y min the two solutions X± coalesce to the value in Eq. (A7). 



In conclusion we have established the equations that describe the position of the singu- 



larities brought by the ^-functions in (21) in terms of the variables X or Y . 



Appendix B: Study of the solutions of Eq. (27) 

With the notations 



u q2 = cos 9 q2 , y = cos6 N , 5<f> = <f> q2 - <f> N , (Bl) 



Eq. (27) for cos 8<p becomes 



Y = u q2 y < vl - y 2 cos 50 . (B2) 

The parameters u q2 ,y,Y are subject to the conditions \Y\, \y\, \u q2 \ < 1. The solution 

54> = ±arccosp(u 92 ), p(u q2 ) = Y Uq ' 2 f = =. (B3) 

v 1 -<v 1 -y 2 

is acceptable, if p has the modulus less than unit. From the expression of its derivative 



dp Yu q2 - y 



(B4) 



du q2 (i_ U 2 2)3/2v /r^2 

we see that for \Y\ < \y\, p is a monotonic function, taking values between —1 and 1 when 



u q2 takes values in the interval 



u q2 e [ u q2Mn = yY- y/l - yVl - Y\ n 92 , max = yY + ^/l- y Wl-Y 2 ] . (B5) 



In terms of angles this condition becomes Eq. (30). If \Y\ > \y\ , then p{u q2 ) has an 



/ y 2 2 

extremum equal to sgn(F)w x ~\ for x = Y/y; the condition \p{u q2 )\ < 1 leads to the 



same interval (B5) 



When equation (27) is solved for the unknown u q2 at fixed (p q2 , we obtain: For \Y\ > \y 



there are two solutions 

(± ) Y y ±^T^\ C os5^yi + {l- y z)pZ-Y* 2 m . 

= + (1 — y ) cos d(f> (Bo) 



y 
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acceptable only for cos 50 E [ J Y ^_^ , 1 ] if Y > 0, and for cos 50 E [— 1, — - 



i-ir 4 



if 



K < 0. For |y| < |y| there is only one solution for any 
or u& if y cos 50 < 0. 



G -7T,vr : m, 



(-) 

92 ) 



if ?/ cos 50 > 0, 



Now we can write explicitly the domain T>(Cn) introduced in Section III A, defined 



by the condition cosa^ < V&n- We refer to the solutions (B6), assuming that in their 



expression Y was replaced by y/C~N > 0. Then the domain V{Cn) can be described by: 

cos 9 g2 E I //. 



'92 



!>0, 07V + 00 



(-) „,(+) 



92 ) U Q2 J J 



where 



arccos 



Cm— COS 2 ( 



sin 2 9 N ' 



for 



cos^at| < a/Cat; 



W — 7T, 0j\r + 7r] , COS 6q 2 E [Uq 2 \ 1] if COS 



'92 



COS 



y 92 



!>/\r) > 0, for cos^/v > \/Cjv; 



92 



— 7T) 0jV t- •'. ] , wo <^ 2 c L — x ' "92 



COS 



"92 



+ 7r], COS o „ G [— 1, M 

> for cos9 N < — v / C ? a7- 



(-)l 



if cosf 



J 92 



"92 



> and cos ^ G if 



> and cos 6 q2 E [—l,u& ] if 
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